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Due to increasing research being conducted at ICASE in the field of fluid mechanics, 
future ICASE reports in this area of research will be printed with a green cover. Applied 
and numerical mathematics reports will have the familiar blue cover, while computer science 
reports will have yellow covers. In all other aspects the reports will remain the same; in 
particular, they will continue to be submitted to the appropriate journals or conferences for 
formal publication. 
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ABSTRACT 

The weakly nonlinear interaction of oblique Tollmien-Schlichting waves and longitudinal 
vortices in compressible, high Reynolds number, boundary-layer flow over a flat plate is con- 
sidered for all ranges of the Mach number. The interaction equations comprise of equations 
for the vortex which is indirectly forced by the waves via a boundary condition, whereas 
a vortex term appears in the amplitude equation for the wave pressure. The downstream 
solution properties of interaction equations are found to depend on the sign of an interaction 
coefficient. Compressibility is found to have a significant effect on the interaction properties; 
principally through its impact on the waves and their governing mechanism, the triple-deck 
structure. It is found that, in general, the flow quantities will grow slowly with increasing 
downstream co-ordinate; i.e. in general, solutions do not terminate in abrupt, finite-distance 
‘break-ups’. 


'This research was supported by the National Aeronautics and Space Administration under NASA Con- 
tract No. NAS1-19480 while the author was in residence at the Institute for Computer Applications in 
Science and Engineering (ICASE), NASA Langley Research Center, Hampton, VA 23681. 
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1. INTRODUCTION 


The nonlinear interaction between two oblique three-dimensional Tollmien-Schlichting 
(TS) waves and their induced streamwise (longitudinal)- vortex flow is considered theoret- 
ically for a compressible boundary-layer flow; this study is an extension of the paper by 
Hall & Smith (1989) who considered an incompressible boundary-layer flow. The same 
theory applies to destabilisation of an incident vortex motion by sub-harmonic TS waves, 
followed by interaction. The interaction is considered for all ranges of the Mach number 
in order to investigate the effect of flow-compressibility. 

The motivation for such a study is essentially the same as expressed by Hall & Smith 
in the introduction to their paper; namely that often in experimental studies of laminar- 
to-turbulent transition on a flat plate (eg. Aihara & Koyama, 1981; Aihara et al , 1985), 
there appear to be longitudinal vortices co-existing, and interacting, with the viscous TS 
modes. As there is no concave curvature of the surface, these longitudinal vortices are 
not Taylor-Gortler vortices driven by surface-curvature (see Hall, 1982a, b and subsequent 
studies); instead one could postulate that they are in fact being driven by, and/or inter- 
acting with TS modes. The reader is referred to the paper by Hall Sz Smith (1989) for a 
fuller account of relevant experimental findings, as well as supporting computational work 
(see, for example, Spalart &: Yang, 1986). These experimental studies are all for incom- 
pressible flow; the author is unaware of any experimented work specifically relevant to this 
compressible study. We note in passing that, for compressible flow over a heated plate, 
buoyancy-driven longitudinal vortices may also be possible (see Hall Sz Morris, 1992). 

Recently, the origin of streamwise vortices in a turbulent boundary layer has been 
investigated theoretically by Jang et al (1986). The Reynolds number is taken to be 
finite and their formulation is of the Orr-Sommerfeld-type. They show that two oblique 
travelling waves can combine nonlinearly to produce a stationary, streamwise vortex — 
this is essentially the theoretical idea later used by Hall & Smith in their independent work. 
However the latter’s approach, that adopted in this paper, takes advantage of the feature 
that the Reynolds numbers of interest in reality are large and so the Reynolds number 
is taken as a large parameter throughout. The nonlinear interaction is powerful, starting 
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at quite low amplitudes with a triple-deck structure for the TS waves but a large-scale 
structure for the induced vortex, after which strong nonlinear amplification can occur. 
Non-parallelism is accommodated within the scales involved. 

The nonlinear interaction is governed by a partial-differential system for the vortex 
flow coupled with an ordinary-differential equation for the TS-waves’ pressure. The solu- 
tion of this coupled system depends crucially upon so-called interaction coefficients which 
are functions of the Mach number; additionally, the TS waves are significantly affected by 
the inclusion of compressibility. It is found that the interaction coefficients, for subsonic 
flow, do not differ significantly in nature from the incompressible ones, but as the flow 
becomes supersonic the restriction (for high Reynolds numbers) that the TS waves must 
be directed outside the local Mach-wave cone (Ryzhov, 1984; Zhuk h Ryzhov, 1981) ex- 
cludes a particular flow solution which is only possible for less oblique modes. The flow 
properties point to the second stages of interaction associated with higher amplitudes. 

It is found that the present formulation breaks down as the Mach number becomes 
large: for then, even when the presence of shock/boundary layer interaction is neglected, 
the viscous sublayers coalesce to form a single boundary-layer. The structure applicable 
in this hypersonic limit has been considered by Smith (1989) and Blackaby (1991). 

The theoretical idea is basically that, if two low-amplitude TS waves are present, 
proportional to E\ t 2 = exp[i(«X=bf Z— ftf)] say; then nonlinear inertial effects produce the 
combination E 1 E^ 1 = exp [zb] , = £3 say, at second order, among other contributions, i.e. 
a standing-wave or longitudinal- vortex flow is induced. Here a, /? and Q are real-valued 
scaled wavenumbers and frequency, whilst X, Z and t are scaled length- and time-scales 
(see later). Equally, the combination of the vortex and one TS wave provokes the other 
TS wave. 

Since the Reynolds number is assumed to be large, the TS waves are supported by the 
triple-deck structure (Smith, 1979,89); however an extra sub-boundary layer and a further 
streamwise length-scale are necessary to capture their interaction with the longitudinal 
vortices. The present vortex/wave interaction mechanism is very similar to that of Hall 
& Smith (1989); the difference is caused by an error in the latter, uncovered by Smith 
&; Blennerhassett (1992). The amended interaction still has the induced vortices lying 
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at the top of the lower deck but now the forcing from the TS waves is solely from an 
inner boundary condition. The wall-shear of the induced vortices modifies the wall-shear 
of the basic flow at the same order as the latter’s leading-order nonparallel correction. 
These corrections to the wall shear force secondary TS waves in the lower-deck, whilst the 
amplitude of the primary TS waves here is governed by an amplitude equation involving 
these corrections to the wall shear. The behaviour of the primary TS quantities at the 
top of the lower deck then leads to longitudinal-vortex activity being forced there. Thus 
the system is truly interactive: the longitudinal vortices are driven by the TS waves, the 
amplitude of which is determined by an amplitude equation involving a vortex-term. 

We consider the interaction for the case of compressible laminar flow over a semi- 
infinite plate. In the next section the underlying boundary-layer flow is outlined and 
the triple-deck structure, for such compressible flows, is reviewed. In §3 the interaction 
equations are derived, and a few special limiting cases of these equations are considered. 
In §4 numerical results are reported and in the last section some conclusions are drawn. 

2. FORMULATION 

We consider the boundary layer due to high-speed uniform flow of a compressible 
over a flat plate. Suppose that L is the distance from the leading edge, and u^, a^, p*^ 
and p^, are the velocity, speed of sound, density and shear viscosity of the free stream 
flow, then we assume that the Reynolds number, Re = is large. This is not 

unreasonable as one is already assuming the presence of a boundary layer. The second 
important parameter is the Mach number, M ^ which we take to be 0(1) for 

the time being. 

A nondimensionalisation based on coordinates Lx (where x is in the direction of flow 
and y is normal to plate), velocities time Lt/u*^, pressure pj^u^p, density 

temperature T^T and shear viscosity p^p is adopted, where the subscript oo denotes the 
value of the quantity in the free-stream. We assume that the fluid’s viscosity and tem- 
perature are related by Sutherland’s formula. Full details of the Navier-Stokes equations 
equations for compressible flow; the resulting boundary-layer equations and associated 
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similarity solutions can be found in several books and articles (eg. Stewartson, 1964). 
Note that only the ‘wall-values’ of the steady boundary-layer flow solution occur within 
scales considered in this paper; however, these quantities depend on the choice of viscosity- 
temperature relation as well as other factors such as whether the plate is cooled and the 
nature of the external pressure gradient. 


2.1 The 3-D compressible triple-deck equations. 


The underlying structure, of the vortex— wave interaction to be considered later, is that 
of the three-dimensional, compressible Tollmien-Schlichting (TS) waves at large values 
of the Reynolds number, namely the three-dimensional ‘compressible triple-deck’. This 
structure has been studied by, in particular, Zhuk & Ryzhov (1981), Ryzhov (1984) and 
Smith (1989); the two-dimensional, compressible triple-deck theory was first considered 
by Stewartson & Williams (1969). Recently, Cowley & Hall (1990) and Duck & Hall (1990) 
have shown that the theory can be adapted to include the effects of a shock for flow over 
a wedge, and cylindrical geometry, respectively; whilst Seddougui, Bowles & Smith (1991) 
have considered the effects of wall— cooling. For definiteness, we assume that the flow is 
supersonic (Moo > 1) during the formulation of the interaction equations; the subsonic 
and other cases follow very similarly and in §4 results are also presented for these cases. 

In the following scalings, the Reynolds number is assumed to be large whilst the other 
factors are taken be be 0(1). The latter are introduced to normalise the resulting governing 
equations as far as possible; however, the Mach number still remains in the upper— deck s 
pressure-disturbance equation and hence it appears in the TS-eigenrelation. 

The stream wise, spanwise length-scales and the time-scale, for Mx, > 1, are 

(x-x 0 ,z-z 0 ) = Re~$Ki(X ,Z ), t = Re”? tJ(M^ - l)"k 

K , = A^/^ritM^-ir*; (2.1a -c) 


here (x 0 , z 0 ) corresponds to the location of the initial disturbance of the laminar base-flow. 
In the viscous sublayer, or lower deck, 

y = Re-SA^'i/i tI(MI - 
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{u,v,w) = Re-^xl^lT^iM^-iyhu, Rt~^ V, W), 

and 

p - Poo = Re-hlpi - 1 )-ip. (2.2a - e) 

The resulting lower deck equations are 

Ux+V Y + W z = 0 , 


Ui + UUx + VU Y + WUz = -Px + Uyy , 

Wi + UWx + VW Y + WWz = -Pz + Wyy, 

P y = 0, (2.3a - d) 


to be solved subject to 

U = V = W = 0, on y = 0; U -* (Y + A(X, Z, t)), as Y -» oo. (2.3 e - /) 


The main deck has 

y = i?e"i /ilriy, (2.4) 

and merely transmits the small displacement effect, A, across the boundary layer as well 
as smoothing out an induced spanwise velocity. The displacement, A, is related to the 
pressure, P, via a pressure— displacement law stemming from matching solutions across the 
three decks (see Smith, 1989). 

The upper deck occurs where 

y = Re-^K,(Ml-iyh-, (2-5) 

note that the Mach number can be scaled out of all but the upper deck equations. 


2.2 Linear Tollmien-Schlichting modes. 

The eigenrelation, for linear supersonic TS-modes, is easily derived using the triple- 
deck scales and equations discussed in the previous subsection. It can be written in the 
form 

(ia)I(Q 2 + y) = (iti7«)(to >{ 4(M f - 1) “ (2 ' 6a) 
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where Ai signifies the Airy function, 

f°° .32 

k= / Ai(q)dq and (o = — i^ft/cv*. (2.66, c) 

hu 

Here a, /i/2 and ft are the scaled wavenumbers and frequency of the mode. The vortex- 
wave interaction to be described concerns only neutral modes, so /i and ft are all real. 

In Figure 1 we present the ‘neutral 1 solutions of the eigenrelation (2.6a), and its 
subsonic counterpart, for a few (illustrative) choices of the Mach number Here (and 

hereinafter) the 4 wave [obliqueness] angle 1 is defined by 

0 < 6 = tai r l (0/2a) < 90°. (2.7a) 

We see for subsonic values of the Mach number (M 0 0 < 1) that neutral modes are possible 
for all wave-angles. However, for increasing supersonic Mach number values (Moo > 1) 
the solution properties start to differ noticeably, with only an ever decreasing range of very 
oblique TS-wave propagation angles, being possible. Thus the restriction (2.7a), which 
can be re-written 

6 > tan -1 [(M£> — 1)^] , (2.76) 

is clearly evident in this figure. We shall see later, once the interaction has been formulated 
and numerical values have been calculated for the important interaction coefficients, that 
this restriction proves to be a more significant ‘compressibility-effect’ on the interaction 
than the ‘direct’ effect due to the Mach number appearing in the interaction coefficients. 

Smith (1989) gives a comprehensive account of the consequences of the eigenrelation 
(2.6a) on the stability of the flow to linear TS-modes (note the factor of 2 difference between 
the definition of /? here and that used by Smith, 1989 and Blackaby, 1991) Our concern 
in this paper is with a vortex-wave interaction based on these length- and time-scales. In 
the next subsection we deduce the size of additional x- and y- scales necessary to capture 
this interaction. 


6 



2.3 The interaction scales. 


In deriving the interaction scales, the same argument as Hall & Smith (1989) is fol- 
lowed but based on the compressible triple-deck scales quoted in §2.1. However, we take 
the coupled lower- and upper-deck equations as our starting point, rather than returning 
to the compressible Navier-Stokes equations; this approach appears simpler. 

We know that TS-waves are governed by the triple-deck structure, and in particular 
by the unsteady interactive boundary-layer equations holding in the lower-deck coupled 
with the upper-deck equations via a pressure-displacement law. If the three-dimensional 
(3D) TS-wave amplitudes are comparatively small, say of order /i « 1 relative to fully 
nonlinear sizes, then nonlinear inertial effects force a vortex motion at relative order /i 2 ; 
the TS-modes are taken to be proportional to 

E x = exp [i(aX + ^Z- (It)}, E 2 = exp[*'(<*X - - (It)} (2.8a, b) 

and we see that combinations yield, in particular, induced longitudinal-vortex terms pro- 
portional to 

E z = exp [i(fiZ)\, (2.8c) 

having only spanwise dependence. 

It can be easily shown that spanwise inertial effects (such as the ‘17 Wx' term of the z— 
momentum equation) decay slowly like 1/Y 2 resulting in the spanwise velocity component 
of the induced vortex to grow logarithmically like In Y (Hall & Smith, 1984,89). Hall 
& Smith (1989) introduced the concept of a new sub-layer (‘the buffer layer’) situated 
within, and at the top of, the lower-deck, along with a longer length-scale (for amplitude 
modulation) to dampen down this logarithmic growth. They showed that the main vortex 
activity was confined to this region. 

Before deriving sizes for the modulation length— scale and the thickness of the buffer- 
layer, we briefly mention the link between the x — scales present and nonparallel effects. The 
triple-deck is a local structure located at nondimensionalised distance x = x 0 from the 
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leading edge. It is short, its length being 0[Re % K\) compared to the 0(1) development 
of the underlying boundary layer, and all the X — dependence of the TS-modes is taken to 
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be in the E\ and E 2 factors. The modulation of the modes is assumed to be on a longer 
x— scale, leaving the eigenrelation (2.6a) unaffected. We define this modulation x— scale, 
X say, by 

x — xq = S 2 X 4~ -Re Re £ K\ S 2 1 y (2.9 u) 


where £2 is to be determined. Thus we have multiple-scales in x; formally we should make 
the replacement 


d d Re t K\ d 

dx~^dx + I 2 W 


(2.9 b) 


in the triple-deck equations (2.3). 

At leading order the wall shear A„, = A^x) is constant (with respect to the X— and 
X— scales) but here we wish to balance the next order term into the interaction equations; 
in fact, for later convenience we have scaled the leading order value, A u ,(x 0 ), out of the 
triple-deck equations. As a Taylor expansion about the local station x = xq, 


A w (x) = A„,(x 0 )(l + ^ 2 XA 6 (xo) H ); (2.10) 


here A/, = A~ l dA w / dx is 0(1) and represents the first influence of nonparallelism (stream- 
wise boundary-layer growth). 

In the buffer layer, where Y = 8Y say (A' 1), the size of the spanwise velocity 

of the induced vortex in the buffer layer is 0(/i 2 In F), ~ h 2 In 8, leading to an induced 
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streamwise velocity of order 


h 2 In 8, by continuity (and noting that the modulation 


X, 


is on X), which alters the basic shear Au, by a relative amount of order 
this is the same order as the ‘non-parallel’ A {,-term if 


8 2 Re^ h 2 

X, T 


In 8 and 


8‘1 


8 2 Re* h 2 


In 8 


(2.11a) 


X{ 8 

The X —modulation has been introduced to damp the induced-vortex velocity components 
in the buffer layer, and this requires the inertial operator, Yd/dX, to balance the viscous 
one, d 2 /dY 2 , i.e. 

( 2 . 116 ) 


3 8 2 Re % 

o <- • 


X! 
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One further relation (between the unknowns /i, 8 and 8 2 ) is required and results 
from balancing the slower X— modulation with the nonparallel effects too (i.e. in the x — 
momentum equation, balancing the A t, term with P~x)i yielding 

S 2 = Re'TSKib . (2.11c) 

The other two sizes now follow immediately from (2.11a-c): 

1 

6~Re'i*Krt and h ~ (Re~T% / \n8y . (2.11 d,e) 

The logarithmic factor in (2.11a,e) is important (Smith & Blennerhassett, 1992); it is 
wrong to dismiss it like Hall Sz Smith (1989). We note that 8 is large, whilst h and 82 are 
small, as required. 

3. THE INTERACTION EQUATIONS 

The method of deriving the interaction equations is identical to that used by Hall & 
Smith (1989); essentially, it involves a standard weakly-nonlinear triple-deck analysis but 
slightly complicated due to (i) the wall shear being weakly z— dependent, and (ii), the extra 
(buffer) layer. Thus here we present only the briefest outline of the derivation; our main 
concern in this paper being the effect of compressibility on their solution properties. Fuller 
details can be found in Blackaby (1991) (see also Hall & Smith, 1984,89; and Smith, 1989). 

3.1 The equations and their derivation. 

The interaction equations comprise of a set of equations for the vortex-terms in the 
buffer-layer coupled with an equation for the TS-wave pressures resulting from matching 
the solutions for the waves between the lower-deck and the upper-deck. The vortex 
equations are forced by the wave via a boundary condition whereas a vortex term explicitly 
enters the wave-equation. 

In the buffer-layer, the vortex terms, u 33 (X ,Y),v 33 (X ,Y) and w 33 (X ,Y), have the 
following sizes relative to the lower-deck scales (2.2) 

(U,V,W) vonex = h 2 In 8 ( Reh 2 u 33 /K u 8v 33 , w 33 )E 3 + c.c.; (3.1) 
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here the notation of Hall & Smith (1989) has been adhered to and c.c.. denotes complex 
conjugate. These scales, together with the wall-shear term SY in the [/-expansion, lead 
to the vortex equations 

U 33 Vyy — J W 33 ~xy ~ *7^^33) 

w 33 yy ~ — 0, (3.2a, b) 

which must be solved subject to the boundary conditions: 

ty 3 3 (X,oo) = u 33 y(X, 00) = u 33 yy(X,0) = 0 and u>3 3 (X,0) = — ifiK\pn\ 2 -, 

(3.2c - /) 

here, 

K = 1 - (/3 2 /4a 2 ) (3.3) 

and pi i is the amplitude of the TS-waves which we choose to be of equal amplitude. Thus, 
we see that the vortex equations are only forced by the TS-waves via a boundary-condition 
which matches the solution in the buffer layer with that found in the lower-deck. 

The desired equation for the pressure amplitude |pi 1 1 of the primary TS-waves can be 
derived by solving the triple-deck equations for the primary and some forced TS-waves; 
after some manipulation we find that it takes the form 

a IpnPOIx + b ^X\pn(X)\ + c « 33F (X,0)|p„(X)| = 0. (3.4) 


This equation was first derived by Hall Smith (1989); however, here the so-called com- 
patibility coefficients a, b and c are functions of the Mach number. The presence of Moo 
in these coefficients is one of the reasons for the solution properties for compressible flow 
differing from those for the incompressible case. In fact, the for supersonic case 


a = 


2r l 7 D6A - 2 

3a 


- iB ~ 5 


2 B B 

^ TT + * 

L 7 3a^ 


b = 


2r l7 .D£o<*A- 5 555 




3a 

(3.5a — c) 


where 


j’i = 


^(£o) 

Ai'^oY 


1'2 = 


«(&) 


5 = 1 + 


£ Q/C (£°) 7 = a 2 + 111 

Ai'(t 9 y 7 4’ 
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(3.6a - /) 


B = 4 (Ml- 1) ~ ° 2 a " d A = 

The coefficients for other flows can be derived very similarly and, in fact, numerical 
results for the subsonic and incompressible cases are presented in the next section. An 
alternative, less physically motivated, derivation of the interaction coefficients is outlined 
by Blackaby (1991) who considers a generalised TS eigenrelation. The quantitative values 
of these coefficients, and their resulting effect on the interaction properties are considered 
later. The (normalised) interaction coefficients 

•Cir = Real(b/a) and c 2r = Real(c/a), (3.7a, b) 


are crucial to the solution properties; especially the large- A" behaviour. 
It is possible to derive a nonlinear ‘integro-differential’ equation, 


4Pnj , 

dX 


C i r ^ b X + 


c 2r fl 2 K Ai' (0) 

r(|)Ai(o) 


l 


X 


\Pn\ 2 mX -*!>)- |p„| = 0, (3.8) 


for the pressure amplitude |pu| from the previous equations. A similar equation has been 
found by Smith & Walton (1989), in their study of vortex-wave interactions. 


3.2 Possible limiting forms for large-A. 

Let us consider analytically the possible flow solutions for large- X. Hall Sz Smith 
(1989) found four such options for their system of equations; however, the necessary 
amendments to their work render one of these options is no longer feasible, namely that 
of exponential growth. Moreover, there is a swap in the signature required for the cru- 
cial quantity Kc 2r for the finite-distance-blow-up and the algebraic-growth-to-infinity 
eventualities to be possible. Thus, the conclusions, drawn later, for the case of zero Mach 
number are quite different from those found in Hall &: Smith (1989). In §4, numerical 
solutions of the interaction equations will be presented and compared with the large- A 
asymptotic predictions that follow. 
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(i) Option I: Finite- distance break-up. 

Hall & Smith (1989) showed that a possible, ultimate behaviour of the nonlinear 
interactive flow, as X increases, was that of an algebraic singularity arising at a finite 
position, say as X — > X Q ; this option is, in fact, still possible for the corrected system 
of interaction equations but with some changes in the details. The similarity forms they 
proposed are appropriate, apart from that for the pressure. As X —> X Q , the behaviours 
for the interaction quantities must have the forms: 

|pii| ~ (X 0 u>33 ~ (Xo-X)-lw(r)), 

u 33 y ~ (Xo — X)~*f(f}), where fj = Y(X 0 — X)~^ . (3.9 a — d) 

When these forms are substituted into the interaction equations, the resulting similarity 
equations can be solved and we deduce that we require the quantity 


Kc 2r < 0, 


(3.10) 


for this option of finite-distance break-up to be a possible large W state of the vortex- 
wave interaction. Note the change of sign necessary for this option to be possible; this 
change is due to the modifications found necessary by Smith &: Blennerhassett (1992). 

The next option that we consider is less ‘catastrophic’, as far as the laminar flow is 
concerned, with the solution continuing to downstream infinity. 


(ii) Option II: Algebraic response at infinity. 

This option is still also possible with the corrected equations. The flow quantities 
must have the following forms 


Ip 


ii 


x*hn)> 


1033 ~ X 3 W ( rj ), u 33 y ~ X r) = YX 


(3.11a - d) 


as X — > oo. 

It is easy to show, from the integro-differential equation (3.8), that we require 


Kc 2t > 0 


(3.12) 
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for this option to be possible; note that this is also a different result than Hall & Smith 
(1989) found. 

The third large-X option proposed by Hall & Smith (see also Smith & Walton, 1989) 
is that of an exponential growth as X — * oo. This option is no longer possible as it relies 
on a forcing term in the 1033 -equation that is not present in the corrected equations. A 
further option, mentioned by Hall Sz Smith, that is still feasible; is that of decoupling due 
to linearisation. Here the TS pressure disturbance |pn | becomes very small, and the vortex 
flow then grows slowly on its own with downstream variable X from its initial upstream 
state. However, this option is ultimately unstable to the TS-waves since the nonparallel- 
growth term, proportional to A* , will dominate the vortex skin friction u 33 y(AT,0). 

3.3 The transonic and hypersonic limits. 


There are three obvious limiting cases to consider for the value of the Mach num- 
ber; below we consider the transonic and hypersonic limits, Moo — *■ 1 and Moc — ♦ °° 
respectively, whilst the incompressible case M 0 0 = 0 is considered in the next section. 

(i) The Mach number tending to unity. 

In his study of the eigenrelation (2.6a), Smith (1989) investigated various limiting 
cases, including those of Moo —* 1 and M 0 o — * oo; Here the ‘transonic, limit’ will be 
considered; without loss of generality, we suppose that the flow is (just) supersonic and 
define 

m = (M 2 QO - 1)*, (3.13) 


so that m is small. Smith (1989) showed that, in this case, the TS-wave quantities behave 
like 


(a,/?, ft) ~ (m"ta*,m~t/r,m"in*) + ---. 


(3.14) 


Substituting these into the formula for the interaction coefficients, we find that 


(c lr ,c 2r )~m ^(c? r ,C2 r ) + 


-'l ri °2r 


0 ( 1 ), 


(3.15) 
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so that the interaction scales X,Y, the vortex disturbances u 33 y,w 33 and the TS- pressure 


amplitude |p, 1 1 also need to be scaled with m: 


(X, Y , u 33 y , u> 33 , fpt 1 1) = (m 8 X 


^ ir?* 

, m 8 y . 


2 7 13 

m 8 , fx* 3 y , 7TI 8 ti> 33 , m 1 6 |p, , I 


*) + ••■. (3.16) 


These scales and the resulting set of equations can be used to check numerical results, for 
the general supersonic case, by providing a ‘transonic’ asymptote. We do not consider the 
transonic limit any further here; Bowles (1990) has considered transonic boundary layer 
transition but the author is unaware of any vortex wave formulations for the transonic 
regime. 

(ii) The large Mach number limit. 

Another limiting case that Smith (1989) went on to investigate was the so-called 
hypersonic limit when Moo > 1; this limiting case leads to some interesting consequences 
for the whole triple— deck structure. Thus it would be most instructive to consider the 
same limit here, as our interaction structure is, of course, dependent to a very great extent 
on the underlying triple-deck scales. 

First, let us recap those results of Smith (1989) which are relevant here, before going on 
to investigate the result of increasingly large Mach number on the interaction-coefficients, 
equations and length-scales. For M 0 0 1 the main features revolve around the small 

regime 

(a, /?, 0) ~ (Moo^d, Moo*/3, MqJQ) + • ■ ■ , (3-17) 

where d,/§ and f i are 0(1). Since o,/3 and fl appear in the interaction coefficients (3.7), 
it is necessary to rescale the coefficients as follows: 

(c lr ,c 2r ) = Moo 2 (ci r ,c 2 r) + •••, where Cl, c 2 ~ 0(1). (3.18a, 6) 


It also is necessary to rescale the quantities appearing in the interaction equations, 


(X,F,fi„ 7 > 33 ,|p„|) = (MU'X.MU'Y.MU'hty.Mll'izz. 

(3.18c — g) 
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A set of interaction equations for this hypersonic-limiting case can be derived (Black- 
aby, 1991). These appear exactly the same as the general interaction equations which 
obscures the fact that the whole multi-layered boundary-layer structure is radically al- 
tered as the Mach number increases. It was shown by Blackaby (1991), based on the theory 

1 _ 3 

of Smith (1989), that as /* Re'S , the triple-deck streamwise length-scale, Re 5 K\X, 
rises to become 0(1) in size; implying that a normal-mode decomposition is no longer 
rational because the nonparallelism of the underlying, growing boundary-layer is now a 
leading-order effect. Further, in this limit it was shown that the lower-deck thickens to 
coalesce with the main-deck. 

This collapse of the underlying compressible-triple-deck structure, as the Mach num- 
ber increases, will obviously occur for the large Mach number behaviour of the vortex-wave 
interaction being considered. However the large Mach number destiny of buffer-layer (in 
particular, its thickness) and the amplitude-modulation scale remain to be established. 
Intuitively, as the buffer-region is ‘sandwiched’ between the lower- and main-decks which 
merge into a single viscous layer in this limit, we would also expect the buffer-region to 
collapse into the same viscous layer. Similarly, as the modulation-scale is ‘sandwiched’ 
between the triple-deck’s streamwise length-scale (which emerges as 0(1) in this Mach 
number limit) and the 0(l)-length-scale of the underlying flow, we would expect that the 
modulation-scale also lengthens to that of the underlying base flow (as Moo /* Re$). We 
now show that these suspicions are correct, by formally considering the large Mach number 
properties of the scales involved. 

Recall that, in the streamwise direction, we have the multiple scales, 

d x — ♦ d x + <^ 2 1 djf -j- Re s A j dx ; 

necessary to capture the vortex-wave interaction. The quantities K\ and 82 are as defined 
by (2.1c) and (2.11c), respectively. In the large Mach number limit, we have seen that 

_! _ 3 

% Moo 4 whilst dx ~ (X ~ Moo 1 , 
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so that the unsealed length-scales, L w and L v say, of the TS-waves and the modulation 
of the induced vortices, respectively, are 

3 3 3 

L w ~ Re~ »K\M^o <C 1 and L v ~ 82 <C 1. 

i 

The Sutherland temperature-viscosity relation leads to, ; 

AT, ~ Mj , (3.19a) 

and so 

2 7 

Lu, ~ Re~^M ^ , /" 0(1), as Moo /* Re$ . 

As far as the amplitude-modulation scale is concerned, we find that : 

L v Re~feK?Ml ~Re~TEM£ / 0(1), as M^ / Re}. (3.19 b) 

Thus, as predicted earlier, this modulation scale does indeed rise to 0(l)-size in this limit 
of the Mach number. | 

Now, let us consider the buffer-region; it lies at the top of the lower-deck, where 
the lower-deck normal-variable Y = SY and 8 is defined by (2. lid). For large Mach 
numbers, we have found, (3.18d), that the buffer-region is characterised by the location 

_ l _ 

where Y = Moo^Y ~ 0(1). Thus the buffer- region lies where 

V ~ 8MI ~ ReTEK^Ml ~ ReT$M~™ 

-» O(ReJS) ~ 0(m 1), as Moo / Re*. (3.20) 

Note that for large Mach number, the lower-deck variable, Y, also scales on Moo! in fact 

i 

Y ~ Ml — hence from this and (3.20) we deduce that the buffer-layer merges with the 
lower-deck, which in turn coalesces with the main-deck. Thus the three sub-boundary- 
layers, present for M 0 c ~ 0(1), have all merged into one single viscous layer. 

Summarising, when M 0 0 —> Re$ the four-layered, short-scaled structure underlying 
the vortex-wave interaction collapses into the two-tiered, long structure found by Smith 
(1989) and considered by Blackaby (1991). 
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4 RESULTS AND DISCUSSION. 

This study was motivated by the desire to find out what changes to the theory, pre- 
dictions and conclusions of the original work by Hall & Smith, are brought about by the 
inclusion of compressibility-effects. However, ironically, the changes brought about by the 
correction of the former turn out to be more significant. For this reason, and for later 
comparison, the new results for incompressible flow will also be presented in §4.2. In the 
following subsection, we show how the interaction equations can be ‘normalised’ so that 
their solution depends merely on initial conditions imposed and the sign of I\C 2 r ■ Numer- 
ical solutions of the normalised interaction equations are presented for both choices of the 
sign of Kc 2 r- 

4.1 The interaction equations renormalised. 

In §3.2 we considered possible limiting-forms, for solutions to the interaction equa- 
tions as A — * oo, and found that the sign of the quantity Kc 2r was crucial in deciding 
whether particular limiting forms were, in fact, possible. This suggests that the interaction 
equations, (3. 2,3.4), can be renormalised. This being desirable, we investigated further and 
found this was, indeed, the case. 

Writing 

X = |c lr A fc |-*A*, Y = |c lr A(,|“^F*, wm = -ipK\cz\~ l W* , 

\pn\ = \c 3 \-iP* and u 33 y = -fi 2 K [c lT -A fr | — 5 |c 3 | _1 r*, (4.1a - e) 

where 

c 3 = -Kc 2r /3 2 /\c lr \ b \?, (4.2) 

leads to the normalised system 

W^.y. - Y*W^. = 0, Ty * y * - = W* 

and P^. + [sgn(ci r Xb)X* — sgn(A C 2 r )T*(X* , 0)] P* = 0, (4.3a — c) 
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which must be solved subject to initial conditions (at X* = 0), together with the boundary 


conditions 


W*(X*, oo) 


t*(X*,oo) = ry. (X*,0) = 0 and W*{X\ 0) = P* 2 {X*). 


(4.3 d-g) 


Thus, the interaction equations (and hence their solutions) are dependent only on the 
initial conditions imposed; sgn(c, r A 6 ) and sgn(ATc 2r ). In all the numerical calculations 
carried out, it was found that c lr > 0, whilst A 6 < 0 for a growing ‘Blasius-similarity- 
variable-type’ boundary layer — appropriate to the present study, if we assume that there 
is no significant wall— cooling or pressure— gradient effects. ^Ve therefore set sgn(c] 7 -A( ) ) 1 

and, apart from consistent initial conditions, the only parameter remaining is sgn(A'c 2r ). 
Thus, with hindsight, it is not surprising that the (predicted) solution properties for, large- 
X, depend crucially on the value of sgn(ATc 2r ). Recall that earlier, in §3.2, we noted the 
following predictions: 


sgn(A”c 2r ) 


> 0 : Algebraic response, as X — ► oo 

< 0 : Finite-distance break-up, as A — * X 0 , 


(4.4) 


for the behaviour of the solutions to the interaction equations. 

To check these predictions, the normalised system, (4.3), was solved numerically; for 
both possible values of sgn(A"c 2r ). and for different (consistent) initial conditions. The 
large-X (X* > 1) properties of the solutions were found to depend solely on sgn(Rc 2r ); 
the initial conditions were found to affect only the initial development of the imposed 
disturbances. The equations were solved by taking ‘central differences’ in Y* and ‘forward 
differences’ in X * (following the method of Hall & Smith, 1989); the appropriate numerical 
checks were performed. 

In Figures 2a, b, we present typical results for both values of sgn(A'c 2r )- In both of 
these computations the system was initialised at X* = —1 (upstream of the neutral TS 


point) using 

P* = P 0 *, W* = Pq 2 (1 + F r * 2 )exp[-y* 2 ], r* = (l-^-+F* 2 )exp[-F* 2 ], (4.5a - c) 


with Pq = 0.1. Note that this initial state, which is consistent with the interaction 
equations plus boundary conditions, corresponds to a ‘mixed’ wave/vortex state. Moreover, 
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we see from the ‘forcing’ boundary condition (eg. (4.3g)) that admissible initial states 
cannot consist of just TS waves alone; the longitudinal vortices must initially be present. 
It appears to the author that the initial states used by Hall & Smith (1989) (see their section 
5; particularly figures 2-5) are inconsistent with their system of interaction-equations plus 
boundary conditions; they do not appear to satisfy the boundary conditions. In their 
study of vortex/wave interactions, Smith & Walton (1989) do not comment on the initial 
conditions they choose. 

Returning to Figures 2a, b, we see that these numerical results are in full agreement 
with the theoretical large- AT* predictions, (4.4). Thus, in the following subsections, it 
is sufficient to merely calculate values of sgn( AC^,-) in order to determine the solution 
properties for large-X*; these being of principal interest. 

4.2 The incompressible case (M^ = 0). 

In their study, Hall Sz Smith (1989) considered ‘this’ vortex-wave interaction for in- 
compressible boundary— layer flow. They cleverly deduced the scales and formulated the 
interaction; unfortunately, they made two unrelated errors in their analysis, both of which 
have a significant effect on the results and conclusions. The first of these errors concerning 
the missing logarithmic term in the interaction scales (most kindly pointed out to the 
current author by Dr. P. Blennerhassett and Prof. F.T. Smith), leads to a simpler system 
of interaction-equations, as well as leading to changes in the possible largo A’ states and 
the necessary parameter values for them to be possible. The second, the term °A 3zY 
missing from the left-hand side of Hall & Smith’s (1989) equation (3.9b), was spotted by 
the current author and leads to a corrected form for c, and hence, a corrected value for 
the crucial quantity C 2 r - 

The interaction coefficients a, b, are as given by Hall & Smith (1989); whilst the 
corrected form for c was given by Blackaby (1991) (see also the Appendix A of Smith & 
Blennerhassett, 1992). In Figure 3 the new numerical values for the important interaction 
quantities, cj r and C 2 r are plotted, versus TS-wave obliqueness angle 6] recall that, for 
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incompressible flow, all such wave angles are possible. Note that c Jr > 0 for all 9; whilst 
C 2 r has one zero, at 9 ~ 32.21°. Recalling the definition of K, 


K 


= l-(/? 2 /4« 2 ) = l 



> 0 
< 0 


: if 9 < 45° 

: otherwise, 


we see that, when Moo = 0, 


, f-1 : if 32.21° <9< 45° 

sgn( ( ' 2 r) — \ +i ; otherwise. 

Thus, from this last result and the numerical calculations described in §4.1, we deduce 
the following: (i) if 32.21° < 9 < 45° then the solution to the interaction equations will 
‘blow-up’ in a finite-distance; otherwise (ii) the solutions will grow slowly (far slower than 
the linear TS-solutions if there were no vortices present), with amplitudes proportional to 
algebraic powers of X * , as X* -»• oo. Note that these conclusions are quite different from 
those of Hall & Smith (who concluded that the ‘finite-distance break-up’ option was most 
likely, apart from the small range 45° < 9 <~ 50° where an ‘exponential-growth’ option 
was favoured). Thus, the theoretically-exciting ‘finite-distance break-up’ option is now 
the exception, rather than the rule. 


4.3 The subsonic, supersonic and hypersonic cases. 

For subsonic (Mo 0 < 1) and some supersonic (1 < M^ <~ 1.15) flows, the properties 
of the interaction-coefficients were remarkably similar to those found for the incompress- 
ible case i.e. graphs of ci r ,c 2r against 9 appear very similar to Figure 3. However, the 
TS-wave angle restriction (2.7b) is found to have a far more significant effect for ‘more’ 
supersonic flows — essentially it can be regarded as preventing wave-angles that would 
allow sgn(-Kc 2 r) < 0, corresponding to the finite-distance break-up option. This is illus- 
trated more clearly in Figure 4 where the results are summarised; we see that the 9 — Moo 
plane splits into four regions (labelled I - IV, as shown). Region IV corresponds to the 
‘barred’ area, where no neutral TS-modes are possible. We see how the border of this 
region acts as an ‘abrupt cut-off’ to the larger-Moo extent of Region II ( finite-distance 
break-up option). This is so much so that, for Mach numbers above y/2, the possibility 
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of finite-distance break-up has gone. Thus summarising, in the subsonic case the results 
are almost identical to the incompressible case; whereas, in general, the finite -distance 
break-up eventuality is not possible for supersonic flows, mainly due to the severe cut-off 
restriction for large Reynolds numbers. To illustrate the last point, in Figure 5 we have 
plotted C 2 r versus /? and 6 for M oo = 3 — note (i) that there is no zero for C 2 r , and (ii), 
the very oblique wave-angles encountered (so that K is always negative and, hence, Kc-ir 
is always positive). 

The last set of results that we present are for hypersonic flow over a wedge, as consid- 
ered by Cowley & Hall (1990), in which a shock is fitted into the upper-deck (at y = y s , 
where y is the normal-variable of the upper-deck), leading to a modified form of Smith’s 
hypersonic TS-eigenrelation (the reader is referred to the paper by Cowley & Hall for all 
details of the formulation). In Figure 6, we present results for the first (lowest) neutral- 
curve for the case y s = 1; here c*ch , flcH ~ 0(1) correspond to a,/3 in the notation of 
that paper. It is sufficient to note that, in our notation, 


fl Pch n 

— ~ Moo > 1 , 

Za acH 

and so the waves they consider are, in general, very oblique. Of particular interest here is 
the (small) interval where C 2 r > 0, so that Kc 2 r < 0, corresponding to the finite-distance 
break-up option; this is an effect of the shock. No such interval is found for the ‘higher’ 
neutral curves; this interval appears to be a feature of the ‘lowest’ neutral curves only (for 
each choice of y s ) and corresponds to ‘crossing 1 the ‘divide’ acH = ficH- 

Finally, we report that for the ‘hypersonic and transonic 1 limiting cases mentioned in 
§3.3, the numerical results and the predicted asymptotic behaviours (for the interaction 
coefficients) were in extremely good agreement. 


5 CONCLUSIONS 

Many of the conclusions of Hall & Smith (1989) carry over to the present study 
and so we concentrate on compressibility-related aspects here. In this paper it has 
been demonstrated that, within the triple-deck framework (Re 1), pairs of small- 
amplitude Tollmien-Schlichting waves and longitudinal vortices can interact, leading to 
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mutual growth. We have seen that two possible ‘eventualities’, for the downstream evolu- 
tion of the interaction, exist; one in which the solutions grow relatively slowly as X — ► oo; 
whilst the other terminates at a finite-distance in a ‘break-up’. Further, we have seen that 
the latter is no longer possible, in general, for supersonic flows (Re >1). 

The interaction has been considered for all ranges of the Mach number; corrected 
results for the incompressible case have been presented and the main effect of compress- 
ibility is through its impact on the TS waves via their governing mechanism, the triple-deck 
structure. In the transonic and hypersonic limiting-cases the interaction modulation scale 
X must be rescaled; in the transonic limit this modulation scale shortens, whilst in the 
hypersonic, limit the opposite is true. The investigation of such vortex/wave interactions in 
transonic and hypersonic flows (not their ‘limits’) should prove interesting; note that the 
former flow has been studied by Bowles (1990), whereas the latter flow regime has been 
considered by Blackaby (1991). 

Other effects which could be incorporated into the present theory include pressure- 
gradient effects; wall— cooling effects (see Seddougui, Bowles & Smith, 1991); cylindrical 
geometry (see Duck & Hall, 1990) and spanwise-modulation (cf Smith & Walton, 1989). 
Finally, we note that Hall & Smith (1991) and Walton & Smith (1992) consider the pro- 
perties of ‘strongly nonlinear’ TS-wave/vortex interactions corresponding to larger wave 
amplitudes than those considered in this paper. 
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